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The issue of magnetic confinement in magnetic fusion devices is addressed within a purely magnetic 
approach. Using some Hamiltonian models for the magnetic field lines, the dual impact of low 
magnetic shear is shown in a unified way. Away from resonances, it induces a drastic enhancement of 
magnetic confinement that favors robust internal transport barriers (ITBs) and stochastic transport 
reduction. When low-shear occurs for values of the winding of the magnetic field lines close to 
low-order rationals, the amplitude thresholds of the resonant modes that break internal transport 
barriers by allowing a radial stochastic transport of the magnetic field lines may be quite low. The 
approach can be applied to assess the robustness versus magnetic perturbations of general (almost) 
integrable magnetic steady states, including non-axisymmetric ones such as the important single 
oo : helicity steady states. This analysis puts a constraint on the tolerable mode amplitudes compatible 

with ITBs and may be proposed as a possible explanation of diverse experimental and numerical 
signatures of their collapses. 

PACS numbers: 52.55.Fa,05.45.+a,52.25.Gj 
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I. INTRODUCTION 



In the present hot plasma toroidal devices, it is widely recognized that the essential limitation that prevents 
ignition to be reached is due to an insufficient energy confinement resulting from small-scale turbulence. During the 
last two decades, some decisive progress has however been made to increase the energy confinement time through 
the widespread experimental achievement of internal transport barriers (ITBs) [U, 0. In these regimes, there exists 
O a region of the tokamak plasma core where the local ion and electron thermal diffusivities are substantially reduced, 
nearly down to their neoclassical values Q. This is likely to be a signature of the regularity of the central magnetic 
I ■ surfaces. In this way, the emergence of ITBs results in a sort of virtuous circle for fusion with both high central 
| temperature T and high energy confinement time te increasing the triple product n e TTE, which implies a better 
. fusion yield. The importance of the magnetic shear, through the form of safety factor profiles, in reaching these states 
by quenching the core plasma turbulent transport has been firmly established. Decisive ingredients were given notably 
by the experimental results obtained by Eriksson et al. Q who demonstrated experimentally that the g-profile could 
be used as a single control parameter to obtain ITBs. In these experiments, they could obtain ITBs for non-monotonic 
g-profiles but not with monotonously growing q, while other parameters such as the E x B flow shear were kept the 
same. Additional evidence was later given by Sauter et al. [4] through a dedicated experiment. They controlled 
the (/-profile using inductive current to generate positive and negative current density perturbations in the plasma 
center, with negligible input power. In this way, they were able to demonstrate that the electron confinement could be 
modified significantly solely by perturbing the current density profile. Other experiments in various tokamaks showed 
that having an inner region with low magnetic shear was indeed sufficient to obtain ITBs [U, [E[ . 

In reverse field pinches (RFPs) also, improved performance related to magnetic chaos healing Q has been recently 
achieved experimentally, yet through non-axisymmetric states 043- ^ was observed that, by increasing the plasma 
current, the RFP operated some change of magnetic topology through a spontaneous transition towards a single- 
helicity equilibrium strongly resilient to magnetic chaos. 

These various experimental results emphasize the importance of the magnetic structure as they show that some 
configurations - reversed-shear, or low-shear, g-profiles in tokamaks or some single-helicity equilibria in RFPs - enable 
a substantial improvement of the confinement. This is an invitation to consider a purely magnetic approach of the 
problem [l(| that will be pursued here. This amounts to address the issue of the plasma confinement in terms of the 
magnetic confinement, namely in terms of the confinement of the magnetic field lines. The ensuing approach is thus 
obviously a simplified approximation of a more complex, self-consistent, reality involving a larger number of effective 



'Electronic address: marie-christine.firpo@lpp.polytechnique.fr 



2 



degrees of freedom, including in particular the electric field. It is a sort of zero order approach that relies on the 
essential principle of magnetic confinement devices which is to use the magnetic field to confine the charged particles 
composing the fusion fuel in a plasma state (Tlj . In this sense, ensuring the confinement of the magnetic field lines may 
be viewed as a prerequisite to particle and energy confinement. As it is well known, this purely magnetic approach 
may be addressed within an Hamiltonian formalism that results from the universal divergence- free nature of B [l2| . 
An intermediate, first-order, step towards a fully realistic description of the physics of magnetic confinement devices 
would be to treat the test particle motion [l3[ instead of the magnetic field lines. This is left for a future work. 

In the absence of resonances, that will be specified below, the dramatic improvement on magnetic confinement 
induced by a low magnetic shear condition will be recalled in Section [TT] and its link to ITBs and core stochastic 
transport reduction discussed. The improved confinement observed in RFPs in single-helicity states will be shown 
to result from the same mechanism when moving to appropriate action-angle variables. However this optimistic 
scenario will be tempered since, in actual tokamaks, ITBs may eventually collapse. In Section IlIIl it will be shown, 
using some Hamiltonian models for the magnetic field lines, that this low shear condition happens to be strongly 
deleterious to magnetic confinement when approaching resonances by enabling stochasticity for much lower magnetic 
perturbation amplitudes. This will be further analyzed in Section IIVI A discussion of this phenomenon and of its 
possible consequences will be finally given in Section El 



It is well known that the equations of magnetic field lines can be cast in Hamiltonian form. Actually, the field lines 
of a magnetic, B(x), or any other divergence- free field are the trajectories of a Hamiltonian (See e.g. the review by 
Boozer [ij] or Abdullaev's book [II]). 

In tokamaks, a dependance on the toroidal angle breaking the rotational invariance is analogous to the usual time- 
dependence in Hamiltonian dynamical systems. Generically one is thus left with a one-and-a-half degrees of freedom 
Hamiltonian, that can be decomposed into an axisymmetric, integrable part and a perturbation term that may be 
Fourier decomposed in the poloidal and toroidal angles. For tokamaks, the Hamiltonian reads 



with Ho{ip) — j W(ij)')chj)' and W{ip) = l/q(ip) the winding profile, inverse of the safety factor profile. In the 
tokamak configuration, H is the poloidal magnetic flux and ip the toroidal magnetic flux conjugated to the poloidal 
angle 9. 



In real magnetic fusion devices, the magnetic perturbations given by the wave amplitudes H mn (tj)) may be induced 
by a variety of phenomena. They may result from external effects, such as the ripple, or may result from intrinsic 
phenomena such as MHD activity. In order to address the magnetic confinement problem, one may turn to simplified 
models for the form of the perturbed magnetic field. Restricting to the single poloidal number m = 1, taking equal 
phases Xin and some special form of the radial perturbation Hx n {ij)) identical for all n's yields the continuous version 
of the tokamap [l4| . Then, considering the limit of a large number of toroidal modes makes the associated Poincare 
(toka)map appear naturally through the use of periodic S distributions. 

Using the symmetric tokamap framework, it was recently shown in Ref. [10( that the existence of a low magnetic 
shear region produced a drastic enhancement of the confinement of the magnetic field lines. More specifically, the 
existence of some ips such that Hq(i/) s ) = 0, that amounts to W'(ip s ) = 0, was shown to improve significantly the 
robustness w.r.t. magnetic perturbations of the regular core magnetic field lines inside the ip — ip s surface, acting then 
as a transport barrier, in agreement with previous results obtained in so-called non-KAM systems including mainly 
maps with reversed shear profiles (See e.g. [l5l - l20j ). Interestingly, Ref. [Io[ gave evidence that this effect does not 
necessitate q to be reversed shear: the g-profile can also be monotonous with an inflexion point, which agrees with 
the diversity of ITBs experimental scenarios. The beneficial effect of the radial extent of the low shear zone for safety 
profiles with q > 1 was demonstrated. Recently, these results were more rigorously grounded: Theorems of existence 
of invariant circles for non-twist maps were formulated for the symmetric tokamap [2l| using a theorem about small 
shear systems formulated by Ortega [22j . 



II. BENEFITS OF A LOW MAGNETIC SHEAR 



A. Hamiltonian formulation for the magnetic field lines 




(1) 



B. Some results on the benefits of low shear 
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C. Extension to non-axisymmetric states 



It is important to note that the domain of applicability of this favorable regime may be extended from axisymmetric 
equilibria to non-axisymmetric integrable states provided ip is replaced by some suitable action variable. Let us 
consider the case where some magnetic configuration involves not only an axisymmetric component H Q (ip) but also 
some (mo,rio) component whose amplitude A is much larger than that of the other (m,n) waves of order e. The 
Hamiltonian reads then 

H(ip, 9, 4>) = H Q {i>) + AH mano (ip) cos(m - n <p), (2) 

up to 0(s) terms. The dominant part of the Hamiltonian, given in Eq. may be put into an explicitly integrable 
form 

9) = H (mo*) - + AH mono (m V) cos 9 (3) 

through a canonical change of variables, with the generating function i^^, 8, 4 1 ) = ( m o9 — ?io</>)\&, that amounts to 
moving to the (mo,fio) wave frame. Even if H[l{ip) does not vanish, it is then easy to recognize that, provided A is 
large enough, a shearless condition in the action variable 

j = _L I ydO (4) 

associated to the integrable part of the Hamiltonian @ may however be obtained. In this way, there can be non- 
axisymmetric magnetic states that are much more resilient to magnetic chaos than what would be their axisymmetric 
components. An illustration of this phenomenon was given in Ref. Q. It was shown that the magnetic structure 
involving some small m = 1 magnetic island with a separatrix was less robust towards magnetic perturbations than 
a reconnected structure with a higher m = 1 amplitude. This transition in terms of the winding profile in the action 
variable J may be interpreted as follows: In the low amplitude case where a magnetic island is present, the separatrix 
is associated to a vanishing of W{J) within a cusp singularity, whereas, as the amplitude grows, the reconnection 
process eventually makes the separatrix disappears inducing locally the growth of W( J) and smoothing the singularity, 
previously associated to the separatrix, which causes the appearance of a local minimum of W( J). This later favorable 
process may be responsible of the unexpectedly good confinement of the non-axisymmetric self-organized states very 
recently achieved in a RFP Q and in a compact toroidal [23j j plasmas. For the tokamak setting, this non-axisymmetric 
formulation of the low shear or zero- twist condition may open new perspectives to improve magnetic confinement. 



III. DETRIMENTAL EFFECT OF LOW SHEAR IN THE VICINITY OF RESONANCES USING SOME 

HAMILTONIAN MODELS 

Let us now go back to the notations introduced in Eq. ([1}. The above analysis did not take into account the 
possibility of resonances. In particular, the fact that the winding number W s = W(ip s ), associated to the action ip s 
such that Hq(iP s ) — 0, could be in the vicinity of some phase velocity n/m corres pon din g to the (m,n) magnetic 
perturbation was not considered. For instance, works using the tokamap framework [Tol . Ibll |24{ ignored de facto this 
possibility by restricting to the m = 1 case while taking typically g-profiles above 1 (with the noticeable exception of 
the sawtooth- related studies of Ref. [HI). 

In order to illustrate the interplay between low shear and resonances, let us first consider a simplification of the 
system (TTJ) similar to the continuous tokamap case but where only the m = 4 poloidal modes are retained, with 

t dx ~ M 

H(M<I>)= / -r^+KHty) V cos (40 - n4>) . (5) 

For M > 1, this represents a one- and- a- half degrees of freedom Hamiltonian in which the perturbation amounts to 
a superposition of plane waves with winding numbers equal to n/4. It is convenient, and physically relevant, to 
consider the large M case for which the derivation of the symplectic mapping associated to ([5]) is immediate through 
the use of the periodic-5 distribution. For computational purposes, we used a symmetric version [Til . [26| of the later 
mapping that yields a very accurate modeling of the continuous dynamics ([5]) with H (tp) = 2^/(1 + ip). The physical 
motivation to consider this specific form of magnetic perturbation was discussed in Ref. |14j |. Basically, this form 
roughly mimics the radial (or ip) form of MHD tearing modes. It gives a vanishing perturbation in the limit ip — > 
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and ensures that the polar axis cannot be crossed so that, under the map, ip remains positive as it should be. It would 
have been possible to use fully realistic MHD profiles for each (m, n) mode obtained from experiments (as e.g. in the 
study of Ref. [271] ). Our point in the present Section is however to use a simpler framework in order to emphasize the 
specific role of W 3 . 

In order to use W s as a single control parameter, we considered the family of winding profiles W p (ip) — W s —c(ip—ijj s ) 2 
with fixed c and ip s - These correspond to reversed shear g-profiles with zero shear at tp s . To be specific, W s was 
allowed to vary while putting c = .5 and ip s — 0.4 in the whole study. 




FIG. 1: The critical value of the stochasticity parameter K, as defined in the text, as a function of the inverse of the minimum 
of the g-profile, W s , for the parabolic family W p (ip) (dashed line) and as a function of the parameter W a in the linear case 
Wi(ip) (plain line). The vertical dashed lines mark the locations of the resonances W a = 1/4, 2/4 and 3/4. 

In order to quantify the impact of the specific value of the minimum of the g-profile, g m i n , (or of its inverse 
W s — l/<7min), wc defined by K c the value of the stochasticity parameter K above which there exists some initial 
condition with ip < 0.1 whose evolution under ([5]) eventually escapes the physically accessible phase space by crossing 
ip = 1. Practically, for each value of W s and a given value of K, up to one hundred initial conditions were chosen at 
random in the range < ip < 0.1 and allowed to evolve under the symmetric map obtained from Eq. ([3]) for up to 5000 
iterations. If some ip happened to cross ip = 1, the process was relaunched for a smaller K or else, for a larger K, in 
a dichotomic way, so as to converge towards the percolation-like threshold K c value. The obtained results have been 
plotted in Fig. [TJ Let us first briefly comment on the meaning of the indicator K c . Using the form of the magnetic 
perturbations as H(ip) = 2tp/(l + ip) ensures that core magnetic field lines remain regular whereas border field lines 
(at ip = 1) are the first to feel chaos as K increases. Having in mind fusion applications, K c has a practical meaning 
in terms of a measure of the confinement of the magnetic field lines, as it gives a threshold above which magnetic 
field lines starting from the core eventually reach the border. It is also an indicator of the global quantitative measure 
of the stochasticity of the magnetic field lines since, in the present model, there cannot be a connection between 
the core and the border of the ^-domain through a regular magnetic field trajectory. However, with respect to the 
magnetic confinement issue, K c may underestimate chaos, in the sense that there may exist non-negligible subsets of 
the physical central phase space domain in which chaotic field lines are chaotic already for values of K below K c . 

It is straightforward to realize that, for the profiles W p (ip) or Wi(ip) = W s — c(ip — ip s ), the dynamics is invariant 
under a translation in W s of some multiple of 1/4. Effectively the curve giving K c as a function of W s is 1/4-periodic 
in W s up to statistical discrepancies. Some part of the curve obtained for the linear strong shear case Wi(ip) has been 
drawn for comparison for about one period in W s . Basically, the effect of low shear happens to be efficient for W s 
between the ra/4 resonances by increasing significantly K c w.r.t. the high shear case. The Poincare plot of Figure 
[5] corresponds to such a favorable situation. This reflects the fact that, in spite of a non-negligible perturbation of 
magnitude K, winding profiles having a low shear region can sustain almost regular core magnetic surfaces. This 
magnetic chaos healing is likely to explain the stochastic transport reduction in ITBs due to the sole effect of the 
q-profile reported in Refs. [3|, Q. 
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FIG. 2: Poincare plot of the system ([5]) with the parabolic reversed-shear profile W p (ip) for K = 0.03 and W„ = 0.645. 



However, in the vicinity of resonances, this tendency reverses: low shear induces abrupt drops of the magnetic 
confinement whereas quite higher threshold values of K are reached for the high-shear linear winding profile. In this 
case, the statement formulated by Rosenbluth, Sagdeev, Taylor and Zaslavski [28| that strong shear is beneficial to 
preserve the regularity of the magnetic surfaces applies. 

For the reversed-shear winding profile W p (ip), low minima of K c occur for values of W s slightly above the linear 
values of the resonances given by n/4. Figure [3] depicts what happens in this case. The abrupt degradation of 
magnetic confinement in the vicinity of the resonances coincides with the emergence of some local stochastic web 
pattern [29U3~D ] that enables a large diffusion in the action tp. Under the perturbation, separatrices of the island 
chains, a characteristic phase-space pattern of the non- twist Hamiltonian systems, are destroyed and replaced with 
stochastic channels of finite width, visible on this Figure. This produces a path to some deleterious radial transport 
already at low values of the amplitude of the resonant mode. The merging of separatrices is the phenomenon through 
which global chaos enters the system and will be more closely examined in Section IIV1 In another physical context, 
Soskin et al. studied a close phenomenon of nonlinear resonance in the minimal case of a single exciting wave [32| . 
Here, the consideration of a large mode number M, that is physically relevant to account for turbulence in toroidal 
fusion devices, should provide an additional effective randomization of motion in the vicinity of the separatrix. 

All other things being equal, we then considered the Hamiltonian system for the magnetic field lines ((SJ obtained 
by just replacing m — 4 by m = 2 and computed its associated symmetric map. Figure |4] displays the obtained 
values of the threshold K c defined for two cases of winding profiles, the parabolic one W p (i/j) and the quadratic one 
W q (ip) = W s — c(ip — ip s ) , together with the curve previously obtained for m = 4 for the reversed-shear parabolic 
profile. This time K c is plotted as a function of the minimum of the safety factor profile q m i n = 1/W S , a parameter 
more commonly used by the tokamak community. In the m = 2 case, the positive role on magnetic confinement 
induced by the existence of a low shear region away from resonances appears in full clarity. The threshold K c attains 
there very large values and the comparison with the m = 4 case is illuminating: the span in W s between two subsequent 
m = 2 resonances is twice larger than the one between two subsequent m = 4 resonances and K c is then allowed to 
grow to larger amplitudes between resonances. The effect of the amount of flatness around <7 m i n is investigated by 
comparing the results obtained for W p to that for the quadratic flattest profile W q in the m = 2 case. Away from 
resonances, it is clear that expanding the low shear zone improves the resilience of the magnetic structure towards 
perturbations as shown in Ref. However, in the vicinity of resonances, the effect reverses. A lower minima of K c 
is obtained for W q and it is reached closer to the resonance q = 2/n (with n = 1 in Fig. 
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FIG. 3: Poincare plot of the system ([5]) with the parabolic reversed-shear profile W p (ip) for K = 0.03 and W s = 0.52. 
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FIG. 4: Critical value of the stochasticity parameter K , as defined in the text, as a function of the minimum of the g-profile, 
g min = 1/Ws- 

IV. LOCAL ANALYSIS 

We wish to interpret the above results within the continuous Hamiltonian system. Let us consider then the 
Hamiltonian 



* M 
H(4>,6,(j>)= I W(x)dx + KH{^) cos (m 6» 

J n=-M 



(6) 
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for some given poloidal mode mo and take W(ip) = W s — c(ip — ip s ) k , with c > and W s = W(ip s ). Let us proceed 
to a canonical change of variable with the generating function 



F 2 (ip',0,<t>) = (mo0-no</>)V>' + ^ 



(7) 



This yields 



dF 2 

89 
8F 2 

dip' 



= ip = moip' + ip s , (8) 
= 9' = m 9 - n o 0, (9) 
and the new Hamiltonian reads, up to a constant, 

M 

H(rJj',8',(l>) = (m W s -n )i>' ~ T^(m ^') k+1 +KH(m^' + i/>„) V cos (0' + n o - n</>) . (10) 

k + 1 * — ' 



-M 



Sufficiently close to the resonance n = hq, 9' + tiqc/) — n<j) is slowly varying, so that one can average over the fast "time" 
<j). This yields the one degree of freedom Hamiltonian 



Htf, &,<!>) = (m W s - n ) $ 



jfe + 1 



(moV 1 ') + (m ^' + ip s ) cos 6»'. 



(11) 



The study of its topology as a function of W s and K should help to clarify the phase space location where the 
stochasticity, induced by the perturbation to (fTTj). should first emerge as K increases from zero. Figure [5] shows the 
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FIG. 5: Phase space of the Hamiltonian (|11[) for some values of W 3 close to no /mo with mo = 4 and no = 2, for c = 0.5 and 
K = 0.01. 

transition from an homoclinic topology with dimerized islands (plot e) to an heteroclinic topology with twin islands 
(plots g and h) as W s increases while being slightly above no /mo = 1/2 . The bifurcation between these two topologies 
[33| occurs at some reconnection threshold W* (see e.g. Ref. j3J,|35|]) that depends on the perturbation amplitude 
K and can be easily computed by equaling the energy values at the x-points obtained for 9 = and 9 = it. 

As it is well known, the nonintegrable perturbation coming from the other modes will gradually replace, as K 
increases, separatrices by stochastic layers. On the basis of Figure 03 one can argue that the largest extent of the 
stochastic domain shall take place in the case where W s is slightly above no /mo, close to the situation corresponding 
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to the plot e, in which case the stochastic layers coming from the two separatrices shall merge for a sufficiently large 
value of the perturbation amplitude. In this way, this analysis gives account of the fact that the smallest values of 
the critical thresholds take place, for the given ^-form of the perturbation, not for W s = no/mo but slightly larger 
values. The Reader is referred to (2||[36| for extra reference studies on the transition to chaos in Hamiltonian systems 
with non-monotonic frequency. 

V. DISCUSSION 

The previous qualitative, yet systematic, study has emphasized the dual impact of low magnetic shear on magnetic 
confinement in toroidal devices. The main achievement of this approach has been to hopefully clarify the dual impact 
of low magnetic shear through a unified picture. In particular, it has been shown using some Hamiltonian models 
for the magnetic field lines that, away from any consideration on dynamics and stability, there exists potentially very 
low thresholds on magnetic perturbations above which safety factor profiles having a low shear region occurring for 
q- values close to q — m/n induce a loss of magnetic confinement. This may now be used both to predict and to 
interpret experimental observations. 

For instance, assume that some reversed-shear (/-profile associated to an ITB regime evolves in such a way that 
q m i n decreases and crosses e.g. q — 2. Then, if there exists some magnetohydrodynamic (MHD) mode with m = 2n, 
typically a (2,1) mode, having a sufficiently large amplitude, one predicts a collapse of this ITB due to an abrupt 
macroscopic degradation of magnetic confinement. This prediction does not need any consideration on the nature of 
the instability that produced this particular (2,1) mode. However, to make a quantitative statement on the amplitude 
threshold inducing the collapse one should have some knowledge on the radial form of the modes involved . 

This prediction does not either preclude the fact that ITBs may be more easy to trigger experimentally for rational 
values of q m [ n for low MHD activity, a fact on which there does not seem to be a general agreement [H anyway. There is 
some subtlety here: For instance, for the original symmetric tokamap, that corresponds to the case where only m = 1 
(non-resonant) perturbations are retained, the most robust magnetic surfaces would actually be obtained for q m i n close 
to 2. It could however explain the fact that the good confinement phase obtained in these cases is always transient 
(37| as various (m, n) modes should eventually be destabilized either linearly or nonlinearly through mode couplings. 
Moreover, as it is well known from the tokamak linear MHD theory, the low m and n MHD modes are typically 
the most dangerous ones, in the sense that they are associated to the largest linear growth rates. Consequently, one 
expects drops of magnetic confinement to take place in the first place for low-order rational values of (/mm, as the 
amplitude thresholds requested for the large-scale degradation of the magnetic confinement should be first attained 
for those low (m, n) modes. 

As illustrated in Fig. 21 the negative impact of low shear should be either benign and transient in the case where 
the q-profile happens to rapidly cross low rational values of g m i n , or particularly deleterious, if the g-profile happens 
to be clamped and very flat about q m i n close to some q = m/n, with a large enough amplitude in the (qn,n) modes. 
These ITBs collapses triggered by stochastization of the magnetic field lines should manifest through flattenings of 
the pressure and temperature profiles in the stochastic region. 

Interestingly enough, these predictions agree with numerous experimental and numerical observations on ITBs 
dynamics. We should first note that experimental evidences indicate that the so-called Edge-Localized modes (ELMs) 
may not be responsible for ITBs collapses since these are observed to happen already before the ELMs occur [38| . 
Additionally, ITB do collapse while the q-profile is reversed-shear and there is a variety of experimental report of core 
confinement collapses in reversed-shear cases when g m j n reaches values close to low order rational values such as 3/2, 
2 or 3 [38l - [40| . The extensive experimental and numerical results of Ref . [4l| consistently show that pressure crashes 
associated to core magnetic stochasticity are strongly correlated with the closeness of g m i n to 2 and with a sufficient 
amplitude in the (2, 1) mode. The proposed picture is also consistent with recent numerical results showing that an 
ITB collapse occurs for g m ; n in the vicinity of 4/3 when turbulence involves a sufficiently large amplitude of the (4, 3) 
mode 42]. Moreover, the above predictions should be most easily evidenced on electron ITBs, since due to their small 
Larmor radius, electrons provide a scan of the magnetic field structure. In this respect, it would be interesting to 
investigate further the connections with the experimental results of Ref. [43| in which the disappearance of electron 
transport barriers are associated with q m j n crossing simple rational values. 

Finally, resonances with m — n associated to the q — 1 surface are by no way special. Sawteeth crashes could 
then be interpreted as ITB collapses associated to q = 1. A complete modeling would however require an enlarged 
framework that is beyond the scope of the present article. 
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